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Abstract 



The QCD /3-function and the anomalous dimensions for the Coulomb po- 
tential and the static potential first differ at three loop order. We evaluate 
the three loop ultrasoft anomalous dimension for the Coulomb potential and 
give the complete three loop running. Using this result, we calculate the lead- 
ing logarithmic Lamb shift for a heavy quark-antiquark bound state, which 



includes all contributions to the binding energies of the form ma^(as Ina^)'^ 
k>0. 
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I. INTRODUCTION 



In this paper, we construct the three-loop anomalous dimension for the Coulomb po- 
tential in non-relativistic QCD (NRQCD) [|l],|2[. The formalism we use was developed in 
Refs. and will be referred to as vNRQCD, an effective theory for heavy non-relativistic 
quark-antiquark pairs. Part of our computation is related to the running of the static po- 
tential 1]^, however effects associated with motion of the quarks do play an important role. 
Our final result for the Coulomb potential differs from the static potential at terms beyond 
those with a single logarithm (i.e. starting at four loops). Combining our Coulomb potential 
running with previous results for the running of the 1/m and potentials allows us 

to compute the next-to-next-to-leading logarithmic (NNLL) corrections in the perturbative 
energy of a heavy QQ bound state, which includes the sum of terms ma^{as Ina^)'^, /c > 0, 
where m is the heavy quark pole mass. This contribution is the QCD analog of the QED 

In a Lamb shift computed by Bethe. In QED, the series m a^{a In a)^ terminates after the 
k = 1 term [§. In QCD, there is an infinite series due to the running QCD coupling as well 
as the presence of non-trivial anomalous dimensions for other QCD operators. The results 
presented here also contribute to the NNLL prediction for the cross section for e"'"e~ — ti 
near threshold [Q. Implications for hb sum rules will be addressed in a future publication. 

The expansion parameter of the effective theory is the quark velocity v. A quark has a 
momentum of order mv and an energy of order mv'^. We assume that m is large enough 
that mv"^ ^ Aqcd and a perturbative description of the bound state as a Coulombic system 
is valid. For a Coulombic bound state, is of order v and contributions suppressed by 
both V and are of the same order. It is useful to distinguish between powers of and v 
when carrying out the matching and when evolving couplings and operators in the effective 
theory, and to only take u ~ for the power counting of bound state matrix elements. In 
the effective theory, the quark-antiquark potentials appear as four-quark operators [0. A 
potential of the form a^/|k|'^ is of order a^f where k is the fermion momentum transfer. 
With this power counting the time-ordered product of a v"- and v'^ potential is of order v"-^^. 
Up to next-to-next-to-leading order (NNLO) the heavy quark potential has contributions 
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The order as/v ~ 1 term in Eq. (|T]) is the Coulomb potential generated at tree-level. The 
next-to-leading order (NLO) term is the one-loop correction to the Coulomb potential. The 
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NNLO terms are the two-loop correction to the Coulomb potential, the one-loop value for the 
l/(m|k|) potential, and the tree-level contribution to the order Ikl^/m^ potential. Writing 
the QQ energy as E = 2m + AE', the terms in Eq. (|T]) generate contributions of the following 
order in the binding energy: 
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In Eqs. (|I]) and the expansion has been performed at the scale /i = m so the cou- 
pling constants are a^ = as{'m). A typical perturbative expansion contains logarithms of /i 
divided by the various physical scales in the bound state. If the logarithms are large, fixed 
order perturbation theory breaks down, and one finds a large residual /i dependence. One 
can minimize the logarithms by setting yU to a value appropriate to the dynamics of the 
non-relativistic system. This is accomplished by summing large logarithms using the renor- 
malization group, and using renormalization group improved perturbation theory. For QQ 
bound states, the large logarithms are logarithms of f ~ a^, and can be summed using the 
velocity renormalization group (VRG) 0. For the binding energy this gives the expansion 
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where the terms are the leading log (LL), next-to-leading log (NLL), and next-to-next-to- 
leading log (NNLL) results respectively. 

In the VRG, one uses a subtraction velocity v that is evolved from 1 to v. This simultane- 
ously lowers the momentum cutoff scale /i5 = mv from m to mv and the energy cutoff scale 

= mu"^ from m to mv"^. The VRG properly accounts for the coupling between energy and 
momentum caused by the equations of motion for the non-relativistic quarks. QED provides 
a highly non-trivial check of the VRG method. In Ref. it was used to correctly repro- 
duce terms in the subleading series of logarithms, including the a^ In^ a corrections to the 
ortho and para-positronium decay rates, the a^ In^ a hyperfine splittings for Hydrogen and 
positronium, and the a* In^ a Lamb shift for Hydrogen. The difference between the VRG, 
which involves the evolution of the momentum and the energy scale in a single step, and a 
conventional two stage renormalization group treatment, m — > mv — > mv'^, was examined 
in Ref. fT^. 

In section |T^ we review the definition of potentials for non-static heavy quarks in the 



effective theory. In section III we compare these potentials with the Wilson loop definition 
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which is appropriate for static quarks. In section |^ we rederive the leading-logarithmic (LL) 
and next-to-leading-logarithmic (NLL) results for the QQ binding energy using the effective 
theory and discuss the two-loop matching for the Coulomb potential using the results in 
Refs. |Tl|JT^ . We also discuss some subtleties in the correspondence between diagrams 
in the static theory and soft diagrams in the effective theory. In section |V] we compute 
the three-loop anomalous dimension for the Wilson coefficient of the Coulomb potential. 
Results for the NNLL energy are given in section followed by conclusions in section |VI1] . 
In Appendix ^ we give some technical details on the structure of divergences in the effective 
theory, and in Appendix |B] we list some functions that appear in the energy at NNLL order. 



II. THE VNRQCD POTENTIALS 



The effective theory vNRQCD has soft gluons with coupling constant as{v)i ultrasoft 
gluons with coupling constant au{v), as well as quark- ant iquark potentials. The potential 
is the momentum dependent coefficient of a four-fermion operator: 



= - ^ (P' p') kp'^ x-p'"^ x-p 



(4) 
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where spin and color indices are suppressed. The coefficient V (p, p') has an expansion in 
powers oi V, V = V(_i) -|- V(o) + V(i) + . . ., where V(_i) = Vc is the Coulomb potential. For 
equal mass fermions 
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where k = p' - p, S = (cri + 0-2) /2, T(k) = {6'^ - /k'^)(T\(Ti. The Wilson coefficients, 
yC^",!) depend on the subtraction velocity z/. In Eq. the color decomposition V = (T^ ® 
T^)V^'^^ + {1® 1)V^^^ has been used and the potential in the color singlet channel is V^^^ = 
yW — (jpyiT)_ (jYhe Casimirs of the adjoint and fundamental representations are denoted by 
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Ca and Cp, respectively.) At LL order the running of the coefficients V^^'^^ was computed 
in Ref . 0] and Vj^l in Ref . mi3| , while the NLL order running of V^. ' was computed 
in Ref. 0. In this work we compute the running of Vi") at NNLL order. This allows the 
computation of the QQ energy spectrum at NNLL order. 

In vNRQCD additional potential-like effects are generated by loops with soft gluons, 
for which the Feynman rules can be found in Refs. [3,0. Matrix elements of soft gluon 
diagrams contribute to the energy beginning at NLO. In contrast, matrix elements with 
ultrasoft gluons start at N^LO. The renormalization group improved energies are obtained 
by computing the anomalous dimensions for these soft interactions and the four fermion 
operators in Eq. (|^). 



III. THE STATIC POTENTIAL VERSUS THE COULOMB POTENTIAL 

Parts of our analysis are related to the study of the static limit of QCD which describes 
heavy quarks in the m ^ oo limit. We therefore briefly review the pertinent results which 
have been derived in this framework. 

In position space the static QCD potential is defined as the expectation value of the 
Wilson loop operator, 

Vstat(r) = lim ;^ln/TrPexp-z5( / Af,dx''\ , (6) 

T^oo i \ JC / 

where C is a rectangle of width T and fixed height r. This potential is independent of 
the mass m of the quarks and depends only on r. In QCD perturbation theory the static 
potential is known at two-loop order . These calculations use static fermion sources 

with propagators which are identical to those in Heavy Quark Effective Theory |]14|. The 



exponentiation of the static potential [|T^,|l^] guarantees that one can avoid dealing with 



graphs which have pinch singularities in momentum space. The analysis of Refs. |[15|J16 
also gives a prescription for the color weight factors for different graphs based on the c-web 
theorem. 



In Ref. [|T^, Appelquist, Dine and Muzinich (ADM) pointed out that at three loops the 



static potential in Eq. (|^) has infrared (IR) divergences from graphs of the form in Fig. |^a,b. 
In the color singlet channel Fig. ^ has color factor CfC\{Ca — '^Cp) while Fig. [l|b is 
proportional to C'jpC\. Taking into account the exponentiation of Vst&t using the c-web 
theorem, the color singlet contribution to V^stat from Fig. |T|(a,b) is simply Fig. ^ with the 
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FIG. 1. Graphs contributing to the three-loop IR divergence of the QCD static potential. 

color factor replaced by C'^Cp, and is IR divergent. ADM showed that this IR divergence 
could be avoided by summing a class of diagrams — those of Fig. |l] with the addition of an 
arbitrary number of gluon rungs. Summing over these diagrams regulates the IR divergence 
by building up Coulombic states for the static quark sources. The summation gives a factor 
exp ([y(^)(r)-r(°)(r)]r) for the propagation of the intermediate color-octet QQ pair, where 
V^'^{r) and V^°^ (r) are the color-singlet and color-octet potentials. The exponential factor 
suppresses long-time propagation of the intermediate color-octet state, and regulates the IR 
divergence by introducing an IR cutoff scale of order [y*^^)(r) — V^°'^{r)] ~ a^/r. 

In Ref. [0] the ADM divergence in the static potential was studied by Brambilla et 
al. using the effective theory pNRQCD [0. They made the important observation that 
along with potential contributions, the definition in Eq. (j^) contains contributions from 
ultrasoft gluons, and the latter are responsible for the ADM IR divergence. They showed 
that the ADM IR divergence in QCD matches with an IR divergence of a pNRQCD graph 
that describes the selfenergy of a quark- antiquark system due to an ultrasoft gluon with 
momenta g'^ ~ Os/t. Therefore, the static potential in pNRQCD can be defined as a 
matching coefficient of a four fermion operator, as in Eq. (P, in an infrared safe manner. 
We will refer to this potential as the soft-static potential. The ultrasoft pNRQCD graph also 
has an ultraviolet divergence. Brambilla et al. computed the coefficient of this divergence 
and extracted a new ln(/i) contribution to the soft-static potential. In Ref. [0 the three-loop 
anomalous dimension was computed for the soft-static potential in this framework. In the 
color singlet channel for scales fi ~ as{r)/r their solution reads 



Ktat(/i,r) = v;ir^^)| 



3/3o ' Vas(/i) 
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where the first term is the two loop static potential derived in Refs. |TT|JT^ . 

For large but finite m the effective theory for QQ bound states has an expansion in 
V. The quark potential in this case differs from that in the static case, and in general one 
cannot obtain the static theory by taking the m 00 limit. To illustrate this consider as 
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FIG. 2. Loop with an insertion of the l/(m|k) and potentials. 

an example the loop graph involving the iteration of a l/(m|k|) and a Coulomb potential 
as shown in Fig. ^ In the effective theory for non-static fermions, the fermion propagators 
give a factor proportional to m. The final result for the graph is independent of m and is the 
same order in v as the Coulomb potential. On the other hand, if we first take the static limit 
m — > oo, then l/(m|k|) — 0, and there is no such graph. This example illustrates the general 
result that the static theory is not obtained as the m ^ oo (or v 0) limit of the non-static 
effective theory. Loop integrals can produce factors of m or 1/v, and thereby cause mixing 
between operators which are of different ordersQ in v. Furthermore, the energy and the 
momentum transfer in Coulombic states are coupled by the quark equations of motion. We 
will see that for this reason the anomalous dimension of the static and Coulomb potentials 
differ at three loops. The three-loop matching for the static and Coulomb potentials can 
also differ. 



IV. THE COULOMB POTENTIAL AT ONE AND TWO LOOPS 

In this work dimensional regularization and the MS scheme will be used. The MS QCD 
running coupling constant will be denoted by and is determined by the solution of 

the renormalization group equation 

/^^^ =/5(a.(Ai)) 

Oo'^lM. Or'^IM 9/^^^M^ 

In a mass-independent subtraction scheme (3q and (3i are scheme- independent. The notation 
aj"^ (/i) will be used to indicate the solution of Eq. with coefficients up to /3n-i (i.e. n 
loop order) kept in the /3-function. 

In the VRG, the soft and ultrasoft subtraction scales fis and fiu are given by fis = 
and /if/ = mz/^. We define the soft and ultrasoft anomalous dimensions 75 and •ju as the 



^However, we stress that if powers of are also counted as powers of v, then operators 
which are higher order in v never mix into lower order operators. 
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derivatives with respect to In /i^ and In/ij/, respectively. The derivative with respect to 
Inz/ gives the total anomalous dimension, 7 = 75 + 27{/. vNRQCD has three independent 
but related coupling constants that are relevant for our calculation: the soft gluon coupling 
05 (z/), the ultrasoft gluon coupling au{i^), and the coefficient of the Coulomb potential 
Vc{i^). The tree-level matching conditions at {fi = m 4^ u = 1) are 

asil) = au{l) = a^{m), Vf )(1) = 47ra,(m) , V«(1)=0, (9) 

and the solutions of the one-loop renormalization group equations for the coupling constants 
in the effective theory are |[T9|J^ 

Vf )(z/) = 47ra W(mz/) , Vi^\u) = 0. (10) 

In deriving the above equations, it has been assumed that any light fermions have masses 
much smaller than mf ^, so that there are no mass thresholds in the renormalization group 
evolution. If there is a mass threshold larger than mv"^ and widely separated from mv and 
m, then it is possible to also include such effects in the effective theory, see Ref. f^. 
At leading order the Hamiltonian for the color singlet QQ system is 

Ho = - + \P-. (11) 

To minimize large logarithms in higher order matrix elements we run 1/ to the bound state 
velocity vi,, which we define as the solution of the equation 

_ acjiy = Vb) _ aW(mt;b) 
Vb — — , (^IzJ 

n n 

where for convenience we have defined 

ae(^) = -^^, (13) 

and n is the principal quantum number. The LL binding energy is then simply the eigenvalue 
of the Schrodinger equation, Ho\ipn,i) = ^E\ipn,i) with the LL solution for the Coulomb 
potential, Vj,'\u) = -CfVP{u) from Eq. (0). Thus, 

AE^^ = [a,{u)f (14a) 
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where in the second hne we have evaluated the energy at the low scale v = Vh- Higher order 
corrections to the energy are all evaluated as perturbative matrix elements with the leading 
order wavef unctions, \il^n,i)- 

Consider how the results in Eq. (|T0|) are extended to higher orders. The graphs for the 
renormalization of the ultrasoft gluon self coupling have the same rules as for QCD, and 
those for the renormalization of the lowest order soft gluon vertex have the same rules as 
for HQET. Since the momenta of soft and ultrasoft gluons are cleanly separated there is no 
mixing of scales, so the anomalous dimension for as is independent of au and vica-versa. 
Thus, one expects that in the MS scheme asi^u) = ^^(mi/) and a^/ (z/) = ^^(mz/^) to all 
orders. 

However, the coefficient of the Coulomb potential can differ from Airaglmu) at higher 
orders. At one-loop, the only order a'^/v graph in the effective theory is the soft diagram |^[| 




7+/5oln(§)+aiJ, (15) 



where Pq = 11/3(7^ — 4Tpn^/3 and ai = 31(7^/9 — 20Tpn£/9 in the MS scheme, and ni is 
the number of light soft quarks. The divergence in Eq. (|T^) is canceled by a counterterm for 
V^^"^-*, causing it to run with anomalous dimension — 2/3o«|(z^). The remaining terms in the 
soft graph are identical to the one-loop soft-static potential calculation and also reproduce 
the set of a^/^^ terms in full QCD with dimensional regularization parameter /i. The one- 
loop matching correction to V^"'"'^''(l) is the difference between the full and effective theory 
diagrams and therefore vanishes at the matching scale = /i5 = m. 

A correspondence between the soft-static potential calculation and soft order 1/v dia- 
grams is expected to persist at higher orders in as well. The Feynman rules for the soft 
vertices are almost identical to the HQET rules used for soft-static potential calculations. 
There are a few notable differences. In the effective theory it is not necessary to use the 
exponentiation theorem [T^JT^ to eliminate diagrams with pinch singularities of the form 



dg° . (16) 

(gO + ze)(-gO + ze) ^ ^ 

These are automatically removed in the construction of the tree level soft vertices because 
the 1/go factors in the soft Feynman rules do not contain ie's, and in evaluating diagrams 



^Note that the soft loop includes soft gluons, soft light quarks, as well as soft ghosts. 
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(a) (b) (c) 

FIG. 3. Compton scattering graphs that contribute to the soft vertex. 





(a) 

FIG. 4. Examples of vertices involving soft gluons. 

these poles are ignored. For example 0, from matching the full theory Compton scattering 
graphs in Fig. ^ one obtains the soft vertex in Fig. ^ which is proportional to [T"^, T^]/qQ, 
where terms proportional to {T^,T^} have canceled. In the soft-static calculations this 
cancellation instead takes place at the level of the box and crossed box graphs, and is 
guaranteed by exponentiation. For the soft-static potential it is known that at higher orders 
there are contributions from the i7r6{qo) terms that originate from the ze's. It was exactly 



this type of contribution that was missed in the two-loop calculation by Peter |TT[, and was 



correctly identified by Schroder [|I^ . In the effective theory these delta function contributions 
belong to the potential regime and soft-static graphs with this type of contribution are 
reproduced by operators such as the one shown in Fig. where a soft gluon scatters from a 
potential. Matching induces these operators to account for the difference between the full and 
effective theory graphs for Compton scattering off two quarks. Thus, the treatment of ie's 
does affect the correspondence between soft-static and soft graphs. The total contribution 
of the graphs in the static theory with k^^ ~ mv gluons is reproduced in the effective theory 
by graphs with soft gluons. 

A real difference between the soft-static and effective theory calculations is the way in 
which counterterms are implemented. The soft-static potential is defined by local HQET-like 
Feynman rules and all UV divergent contributions from soft gluons are absorbed into vertex, 
field, and coupling renormalization. The renormalization of the four point function is taken 
care of by the renormalization of the two and three point functions. Renormalization of the 
vNRQCD diagrams is quite different because potential gluons are not treated as degrees of 
freedom. The effective theory has graphs with soft gluons and in addition the four-quark 
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FIG. 5. Order a^/f diagrams with potential iterations. The x denotes an insertion of the 
p^/SmP relativistic correction to the kinetic term. 




FIG. 6. Examples of order a^/v diagrams with soft vertices. The vertex with a cross denotes 
an insertion of a one-loop counterterm. 



Coulomb potential operator. The overall divergence in soft graphs, such as the one in 
Eq. ([T5|) , are absorbed by Vc, while sub divergences are taken care of by counterterms for the 
soft vertices and lower order Vc counterterms. The sum of unrenormalized soft-static and 
purely soft diagrams in the effective theory agree. So if these were the only considerations, 
then the description of these effects would be basically a matter of convenience. However, 
divergences associated with ultrasoft gluons can only be absorbed into Vc, so at the level 
that these gluons contribute in the effective theory it is necessary to adopt the four quark 
operator description from the start (i.e. just below the scale m). 

To calculate the NLL energy we also need the two-loop anomalous dimension for Vc, 
which is obtained from the renormalization of order a^/v diagrams. At this order there are 
effective theory graphs with iterations of potentials, shown in Fig. |^, and soft effective theory 
diagrams, as in Fig. Graphs with ultrasoft gluons do not contribute at this order. The 
potential diagrams are finite in the ultraviolet and reproduce the Coulombic singularities 
in perturbative QCD. The contributions from the soft diagrams can be determined from 
the soft-static potential calculations. For the color singlet channel, the UV divergences in 
the soft-static two-loop diagrams were calculated in Ref. [^] and the constant terms in 
Refs. [Tl|Jl^. The sum of unrenormalized soft diagrams has the form 



ia|(z/) Ct 



Pi A + 2/?oai 
e2 e 



k2 47r 

The effective theory counterterm graphs give 
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^2 



2^ + Min(d^ 
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47r L ' e ^ / e 
Taking the sum of Eqs. (|T^ and (|T^) we find that up to two-loops the counterterm for the 
color singlet Coulomb potential has the form 
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4:716 



(4vr) 



(19) 



where Pi = 34C^/3 — ACpTpUi — 20CyiTp?7,^/3. The a^/e divergence is proportional to 
so the two-loop anomalous dimension for V^'^-' is determined by the two-loop MS /3-function, 
and the NLL coefficient of the singlet Coulomb potential is Vjf\u) = 47raJ^l (mz/) . 

The energy at NLL order involves including the NLL coefficient for the Coulomb poten- 
tial, vi') (z/) in Eq. ( |14a|) , and calculating the matrix element of the one-loop order 1/v soft 
diagram between Coulombic states [ac(z/) = — V^*-*(z/)/47r] 



-Cf aliu) ( ^ 



1 



,a„.n(§) 



(20) 



itiCf a|(i/)ac(i^) 



ln( 



n V 



+ ?/'(n + / + l)+7ij 



As expected, at the low scale i/ ~ u^, there are no large logarithms in the matrix element. 
Combining Eq. ( ppj ) with Eq. ( |14a] ) gives the energy valid at NLL order. 
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(21) 



In the next section, the three-loop running of the Coulomb potential will be derived. We 
therefore need the two loop matching condition, and so consider the finite parts for the two 
loop graphs. The sum of renormalized soft diagrams in Fig. ^ is 



ialiiy) Cf 
k2 47r 



(A + 2/3oai)ln(^) 



02 



(22) 



where from Ref. [|T^ the sum of constants in Eqs. (|T^ and (|T8|) is 02 = a^2' + '^2 
456.75 — 66.354r2£ + 1.235r;,| for rig light flavors. The matching coefficient for Vc at the scale 
m is given by the difference between the 1/k^ terms in the QQ scattering amplitude in 
the full and effective theories. It is convenient to analyze the two loop result in the full 



theory by using regions in the threshold expansion The soft region exactly reproduces 
the result from the soft graphs. Furthermore, the potential region exactly reproduces the 
results for the potential graphs in Fig. ^ Thus, the matching correction for Vc(l) is also zero 
at two-loops. In general, a non-zero matching correction appears when there is a full theory 
contribution from an off-shell region such as the hard regime or when UV divergences appear 
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FIG. 7. Graphs with ultrasoft gluons which do not contribute to the running of the Coulomb 
potential. The divergences in a)-e) are canceled by graph f) which has an insertion of the corre- 
sponding Vk counterterm(s) denoted by (g). Graph g) is UV finite. 

in the effective theory graphs. [| In the full theory at two loops there are no contributions 
proportional to 1 /k^ from off-shell regions. The soft effective theory graphs are UV divergent, 
however these divergences are in one-to-one correspondence with UV divergences in the full 
or static theory. Finally, the graphs with iterations of potentials are UV finite. 



V. THREE-LOOP RUNNING OF Vc 

To compute the three-loop anomalous dimension for the Coulomb potential we need to 
evaluate the UV divergent graphs in the effective theory that are order a^/f. We begin 
by considering diagrams with an ultrasoft gluon. In Coulomb gauge we have graphs with 
p ■ A/m vertices as well as the coupling of ultrasoft gluons to the Coulomb potential from 
the operator 

9 • y{T) fABC 

^ = ^^\r /^^ k ■ (^A^) , T^^p ^l^, f\^^ . (23) 

All graphs with two p ■ A vertices are UV finite or are canceled by two loop graphs with 
insertions of the one-loop counterterms for V2 and Vr computed in Ref. Q]. The remaining 
diagrams are shown in Figs. ^ and ^. Graphs |^a through |^ have UV subdivergences which 
are exactly canceled by the diagram with Vk counterterms shown in Fig. |^. These graphs 
contain subdivergences that were responsible for the running of the l/(m|k|) potential at 
two-loops . Graph |^ is UV finite. 

The divergent diagrams with an ultrasoft gluon which are not completely canceled by 
a counterterm diagram are shown in Fig. |[ Consider the three-loop graph in Fig. ^ with 



^An example where UV divergences in the effective theory affect the matching is the two- 
loop coefficient for the production current. 
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FIG. 8. Graphs with an ultrasoft gluon which contributes to the three- loop running of the 
Coulomb potential. 



momenta / for the loop with the ultrasoft gluon and k and q for the remaining loops. After 
performing the and g° integrals by contours, the loop integration involving / is 



(24) 



From this expression we see that the ultrasoft momentum I and potential momentum, k 
and q, are not completely separable since they appear in the same propagator. The / 
integration produces an UV divergence while the remaining integrations are UV and IR 
finite^. Evaluating the remaining integrals gives 
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Fig. 8a 
where the color factor is 

(Csa 
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(25) 



(26) 



and for gauge group SU(A^e), Cd = N^-A/N^ and Ci = {N^-1)/{AN^). The graph in Fig. ||c 
involves the iteration of a 1/k^ potential and a Vk counterterm and also has a Coulombic 
infrared divergence. This graph cancels the corresponding product of IR and UV divergences 
arising in Fig. IHb. The sum of graphs in Fig. R]b,c still has an UV divergence, and we find 



4z 
T 



Fig. 8b + |c : 
where the color factor is 
(Csfec) = Ca 



(47r)3 k2 



2 

k2. 



32 



(27) 



(28) 



The sum of divergences in Eqs. (^) and (pTQ are canceled by a three-loop counterterm for 
Vc- Differentiating with respect to In/ig and In/i^/ gives the anomalous dimensions 



^For static quarks this three-loop graph also has an IR divergence []T8[, but in the non-static 
case we find that this divergence is regulated by the quark kinetic energy. 
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Ca{Ca + Cd 



32 



\vp{ 
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75 = y 



2CaCi + 



Ca{C/ 



Cd) 



32 



(47r)3 
Vf)(z 



(47r)= 



(29) 



The total anomalous dimension from the ultrasoft diagrams is 7 = 27^7 + 75, so 



32 



2CaCi 



Ca{Ca + Cd)' 



32 



Vf)(z 



(47r)^ 



(u) (t ® f ) . 



(30) 



The presence of an ultraviolet divergence in the ultrasoft graphs induces an additional ul- 
traviolet divergence in the soft graphs (as discussed further in Appendix 0). This divergence 
induces an additional contribution to the soft anomalous dimension: 



7s 



2CaCi + 



Ca{Ca + Ca 



\2n 



32 



a. 



M {T®f) . 



(31) 



For the Coulomb potential the remaining UV divergences in the soft graphs correspond 
to divergences which are canceled in the static calculation by field, vertex, and coupling 
renormalization. As discussed before, these divergences give a contribution proportional 
to the three-loop MS /3-function, so for the color singlet channel we have the additional 
contribution 



is) 

Is 



2CfP2 



' (47r)2 ■ 

For QCD, P2 = 2857/2 - 5033Ti^/18 + 325n2/54 for light flavors. 

Combining Eqs. (|30D-(|32D in the color singlet channel and using the LL relation V^^-'(z/ 



(32) 



4:TTas{i^) gives the total anomalous dimension for the Coulomb potential to three-loop order 



Poasi^)+/3i- 



(4vr) 



—CaCf 



-al{u)au{i^) - ai;{u) 



Solving this equation with the two-loop boundary condition V^*''(l) 
NNLL result for the running Coulomb potential is 



(33) 

-4:7iCFas(Tn), the 



-ATrCpaf'^ (mu) 



6po 



2z 



where 



aS}^{mv) 
a Am) 



w 



11 

T 

a W(mz/2) 
a]^l(mz/) 



4 



+ 41n(w) 



(34) 



(35) 
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FIG. 9. Graphs with double logarithms that are determined from Fig. ga. 

To compare to the static potential result from Ref. we can expand the running cou- 
plings in as{m). Ignoring the aj^l(mi/) term, the remaining logarithms in the V^^^u) coeffi- 
cient are 

- \ CfC\ at{m) \n{v) + ^ CpCl al{m) - CpC\ a%m) \n%u) + .... (36) 
6 Sir loTT^ 

Expanding the analogous contribution to the static result in Eq. from Ref. [0, taking 
/i = mz/2 and r = 1/ mu, gives 

- \ CfCI atim) Hu) + ^ CpCl al{m) - M CpC\ al{m) \n\v) + .... (37) 
6 47r 7/7r^ 

From Eqs. ( |36D and ( p7| ) we see that the single In(i^) terms agree, but the higher order 
logarithms differ. 

The origin of the difference between the running of the Coulomb and static potentials 
is the relation between scales for the case of moving versus fixed quarks. For the Coulomb 
potential the non-relativistic quarks obey the dispersion relation E = p2/(2m), which re- 
lates the energy and momentum scales, and logarithms of E, p and k cannot be treated 
independently. The anomalous dimension in Eq. (|33|) generates ln^(z/) terms that reproduce 
the double logarithms that appear when diagrams such as the ones in Fig. ^ are evaluated at 
the hard scale m. The graph in Fig. ^ includes a vacuum polarization loop for the ultrasoft 
gluon which only sees the scale mv^, while the soft loop in Fig. [pD sees only the scale mv. 
Even though each of the two logarithms in Fig. comes from a different ratio of low energy 
scales, the modes in the graph must be included at the same time since in the effective 
theory the division between soft and ultrasoft modes (and the multipole expansion) occur 
right at the scale m. Furthermore, for all scales below m the couplings for soft and ultrasoft 
gluons run at different rates. As mentioned in the introduction, this correlation of E, p and 
the momentum transfer k has been tested successfully for bound states in QED [Q. 

In contrast, consider the situation with two static quarks where the distance between 
them, r ~ 1/1^1 5 is held fixed externally, and the energy fluctuations are about E = 0. In 
this case the scales r and E are not correlated. Furthermore, operators with powers of 1/m 
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play no role in the calculation of the anomalous dimension, unlike the Coulombic case with 
an expansion in the velocity. The difference between the static and non-static calculations 
occurs essentially because neither the m — > oo nor the w — limit of the effective theory is 
the same as the static theory. 



VI. THE NNLL ENERGY FOR A QQ BOUND STATE 



The energy at NNLL has contributions from matrix elements of operators of order 

1 , asV, V \ . 



There are contributions from tree level matrix elements, which include the p'^/m^ operator 
(~ f ^), the order v potentials with LL coefficients (~ ttsf ), the order potentials with NLL 
coefficients (~ a^v^), and the correction to the energy from the Coulomb potential with the 
NNLL coefficient (~ <y^/v). We also have the matrix element of the order a^/v one loop 
and two loop soft diagrams in Fig. p. Finally, there are the double insertions of two a'^/v 
soft diagrams (~ af/w^). For simplicity all n£ light quarks are taken to be massless. 
The contributions from tree level matrix elements are [adi^) = —Vjf\i')/4:7r] 



y(o) 



4^3(2/ 



(38) 



^ Xijs5si (1 - Sio) ^is)^j^^ _^ ^{Si2)ijs5si (1 - ^;o) yi;*)^^^^) 



47r/(/ + 1)(2/ + 1) 



47r/(/ + l)(2/ + l) 



/zp!\ 



m 



where the Wilson coefficients V, 
the soft loops are 



An 21 + 1 

is) 



ck2sr At defined in section |T[ The matrix elements of 



+ ^ + fii In 

4 ^ \ \aAv 



+ 21n 



2 ' 



n V 



^{n + / + fl 



lE 



+ iV2(n,/) 



(39) 



TU 



V 



n V 



■mC\a%{v) 



a\ , I . I nv 



4 In 



n V 



, Potti In 
4 V \ar\i^ 



lE 



-2iVi(n,/)+7s 



iVi(n, + yJ + 4iVo(n, /) + 47siVi(n, /) + 7I 
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where the functions No^i^2{n,l), {Si2)ijs and Xijg are obtained from Refs. p2|,^ and are 
summarized in Appendix |B[ The only additional contribution is the result for the NLL 
energy in Eq. ( pT] ) with the Wilson coefficients evaluated at one higher order. 

Again, at the scale u = Vb there are no large logarithms in the matrix elements; the large 
logarithms are summed up into the Wilson coefficients. For v = Vb the energy at NNLL 
order in terms of the pole mass m reads 



LL 



AE 



NLL 



_|_ ^£iNNLL 

mCF(yl{vb)ac{vb) 



(40) 



+ 



m 



m 



4^3(2/ 



2f3o{ip{n + I + I) + IE 
2 



ai 



[advb)] 



3 

An 



+ i^^t"^^^'^)] 12^ 
5io 



V^'\vb) + s{s + l)Vi'\vb) 



21 
(2/ + 1 



An) 



inn{2l + 1) 



Vi'\vb) 



47r/(/ + 1)(2/ + 1) 



V^l'{vb) + 



3(5'i2)ijs5si(l - 5, 



PoN^inJ) 



47r/(/ + 1)(2/ 
2 



Vt\vb) 



^(n + / + l)+7E 



02 



mC]pal{vb) 



ANo{n,l) + A-fEN,{n,l) + -fl 



2Ni{n,l)+-fE 



When this expression is expanded in powers of as at a given renormalization scale the LL, 
NLL, and NNLL predictions for the energy become series in (aslnag) as in Eq. (|^). The 
terms beyond NNLO that are determined unambiguously are those up to maf In a^)^, 
A; > 1. From Eq. (^) we see that up to NLL the series are determined by the running of 
a^(yu) since up to NLL order we have ac{v)/CF = = olg{mv). At NNLL order the 

series is no longer just determined by the QCD /^-function since operators besides the strong 
coupling have non-trivial anomalous dimensions. 

It is well known that the convergence of predictions for AE in terms of the pole mass 
are plagued by the presence of infrared renormalons. If predictions are made in terms of 
a short distance mass such as the MS mass the leading renormalon in the pole mass and 
potentials cancel P^ , P3| and the convergence of the perturbation series is improved. 
A phenomenological analysis, which includes the issue of renormalon cancellation in the 
presence of resummed logarithms, will be carried out elsewhere. 

Past finite order predictions for the QQ energies have typically been made using the 
strong coupling evaluated at the soft scale mv. At LO and NLO this is the natural choice 
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FIG. 10. Comparison of the NNLL binding energy predictions (thick black lines) with the 
NNLO predictions (thin green lines) for the l^Si (solid), 1^5o (dashed), 2"^ 5*1 (dot-dashed), and 
2^ Pi (dotted) states, and for different values of the subtraction velocity v. 

since all logarithms a^(asln/i)'^ and ^^(asln/i)^ in higher order matrix elements are mini- 
mized at /i = |k|. Effectively this choice turns the LO and NLO results into the LL and NLL 
predictions. However, at NNLO this choice of /i may not be optimal since matrix elements 
begin to involve factors of [a^ ln(/i/£')]. The NNLL prediction is not generated by a simple 
replacement rule and instead involves the non-trivial Wilson coefficients Vc,k,2,s,r,K,t{^)- 

In Fig. |TD| we compare the NNLL energy predictions at subtraction velocity v (thick black 
lines) with the NNLO predictions with coupling a^(mi^) (thin green Unes), to illustrate the 
impact of summing the logarithms. The displayed states include 5*1 (solid 

lines), 1 (dashed lines), 2"^ 5*1 (dash-dotted lines), and 2^Pi (dotted lines). As input we 
have chosen mf, = 4.8 GeV for the bottom quark pole mass, a^"^^'*-' (m^) = 0.22, and have 
included three-loop running for the evolution of the strong coupling to lower scales. The 
NNLO and the NNLL energies are equal for = 1, because no logarithms are summed into 
the Wilson coefficients. The summation has the largest impact on the n = 1 S-wave states, 
and in all cases reduces the size of the binding energy. In Table | we summarize for v of order 
a typical quark velocity, v = (0.35, 0.4), values for the NNLO result and the renormalization 
group improved NNLL calculation. Relative to the NNLO results the scale uncertainty in 
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bb state (n^^+^Lj) 


1^5*0 




2^ Si 


2^ Pi 


AE at NNLO (MeV) 


= 0.35 




— Q79 


—449 


—488 




u = 0.4 


-912 


-860 


-389 


-423 


AE at NNLL (MeV) 


= 0.35 


-641 


-614 


-382 


-446 




u = 0.4 


-705 


-678 


-356 


-403 



TABLE I. Comparison of the NNLO and NNLL predictions for the binding energy corrections 
(in MeV) for bb states. The NNLL results include the summation of logarithms not accounted for 
by ctg {mv) . Results are shown for two values of u of order the heavy quark velocity. 



the NNLL predictions is somewhat reduced. 

The summation of NNLL logarithms also reduces the size of the ground state hyperfine 
splitting, but increases the size of the analog of the QED Lamb shift, (the 2^Si — 2^Pi 
splitting). For v = (0.35,0.4) and at NNLO the hyperfine splitting and the 2^Si — 2^Pi 
splitting are 



E{l^Si) - E{l^So) = (68, 52) MeV , 
while at NNLL order we find 

E{l^Si) - E{l^So) = (27, 27) MeV , 



E{2^Si) - E{2^Pi) = (39, 35) MeV , (41) 



E{2^Si) - E{2^Pi) = (64, 47) MeV . (42) 



Summing the logarithms reduces the perturbative contribution to the hyperfine splitting by 
a factor of two. The 2^51- 2'^Pi splitting is fairly sensitive to the value of u. 

In order to examine the size of the logarithmic terms that are summed at NNLL order 
we can expand Eq. ( ^Dj ) in powers of a^. To suppress logarithms proportional to the QCD 
/3-function in the energy at NNLO is it convenient to expand in the coupling = a^(mffc): 



AE 



m a. 



— ma. In a. 



(43) 



Ca 
3 



C2 



_A_ I 4 Ca Cp 2 C 
'Y nf2/ + 1) 



1 

1 n 



3 5/0 C'f 



2n 



{Ca+2Cf) - 



^2 

fi 1 2 ^ F 

maMn a 



^10 Cp 
47r^ I 6n 



n(2/ + l) ' n \2l 
7CaCp6io 6si _ Ca Cp (1— 5^) 

4n/(/ + l) (2/ + 1) 

Ca 



A 



3 n 
13 

2 
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(^4:Xljs+{Sl2)ljs 



Cf] +^{25Ca + 22Cf 



6n 

8n/(/ + l) (2/ + 1) 
+ . . . . 



5Po + 7Ca 



/5o ( 2 Xijs + - {Si2)ijs ) + C'a ( 2 Xijs + (Su) 



Ijs 



The af, a^, terms are not displayed, but agree with the result in Ref. The terms 

proportional to marinas agree with Ref. The ma^ln^Os result is new, as are higher 

terms in the series. Numerically we find 



^ ^ = -0.444 al - 1.595 - 9.73 - 8.56 In a, - 3.41 at In^ a, - 15.5 al In^ a. 

+ ... , 



m 



AE{l^Si 



m 



-0.444 ai - 1.595 < - 8.68 < - 6.80 < In a, - 0.904 a° In^ a, - 12.1 a'. In'' a. 



+ ... 



-0.111 - 0.546 af - 3.07a! - 0.961 aMna, - 0.113 a! In^ a, - 1.73alln^a 



.71^3 , 



m 



s ^^-^ "'s 



+ ... , 



^^^^ = -0.1 11 - 0.644 al - 3.22 - 0.398 In a, - 0.005 In^ a, - 0.752 aj In^ 

-f.... (44) 



We note that for (I^Sq, l^S"!, 2^S'i, 2^Pi) using = 0.35, the Ina^ terms in Eq. (|4D give 
(93%, 85%, 84%, 80%) of the complete sum of logarithmic terms. 



VII. CONCLUSION 

The three-loop anomalous dimension for the Coulomb potential was computed in the 
presence of non-static quarks. In terms of the subtraction velocity u the anomalous dimen- 
sion depends on both ^^(mz/) and ^^(mi/^). Our result differs from the three- loop anomalous 
dimension for the potential for static quarks. This is due to the fact that the energy and 
momentum scales are coupled for non-static quarks. 

The perturbative energies for QQ bound states were computed at NNLL order, including 
all terms of order mQ;^(as In a^)'^, k > 0. The main effect of summing the logarithms is to 
reduce the binding energy for the 1^5*0 and l^Si states by an amount of order a few hundred 
MeV. The effect on states with n > 2 is substantially smaller. 
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APPENDIX A: ULTRAVIOLET AND INFRARED DIVERGENCES 

In this appendix we discuss the structure of UV and IR divergences in the effective theory 
and their relation to the renormalization group running. In particular we explain why it is 
correct to treat all 1/e poles in soft gluon loops as UV divergences. This result follows from 
the fact that we have taken the full propagating gluon field and split it into two parts (soft 
and ultrasoft) which fluctuate on different length scales {mv and mv"^) and have different 
Feynman rules. The information that the two gluons came from a single field is reflected by 
a relation between the IR divergences in soft gluon diagrams and the UV divergences in the 
ultrasoft diagrams. Here we address how this correlation is accounted for in the presence of 
the relation fiu = fi'^/m between ultrasoft and soft subtraction scales. 

For simplicity only diagrams with single 1/e poles and either purely soft or ultrasoft 
gluons will be discussed. Also, real Coulombic IR divergences will be dropped. A general 
diagram with a divergent soft gluon loop then gives a soft amplitude with the divergence 
structure: 

= + — + (7 = + (C - 5) . (Al) 



Strictly speaking pure dimensional regularization does not distinguish between UV and 
IR divergences; however a distinction can always be made either by examining the analytic 
structure before expanding about e = or by also performing the calculation with a different 
IR regulator. The notation in Eq. ( [AID slightly redundant so that C represents diagrams 
involving scaleless integrals (such as tadpole graphs), and A and B represent all other graphs. 
For example, for the Lamb shift in QED for Hydrogen only C is non-zero, while positronium 
also has non-zero A and B. At the same order in the power counting as Eq. ( |Al| ) there is 
an amplitude with a divergent ultrasoft gluon loop which has the form: 

+ (A2) 

^UV ^IR 
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The nontrivial information is that the (C — B) term in Eq. ( |A2|) is determined by the IR 
divergence in the soft amphtude in Eq. ( |A1| ). In general D is independent of C — B since 
ultrasoft graphs are not always proportional to 1/euv — ^/^ir- The infrared divergence in 
Eq. ( [A^ ) will match with an infrared divergence in full QCD. 

The renormalization group running is determined by the UV divergences. Naively, in 
Eq. ( [Al[ ) the e^/y's correspond to the scale m and the e//j's correspond to the scale mv, while 
in Eq. ( |A^ ) the Cf/y's correspond to the scale mv and the ej^s correspond to scales of order 
mv"^. However, examining iA^ + iA^ we see that the {C — B) term in the soft amplitude 
simply pulls up or transports the 1/euv in the ultrasoft graph to the hard scale 0. The 
scale dependence of the coefficient C — B in Eq. (|A2|) does not affect this argument since 
the scale dependence of C — i? in Eq. ( |Al| ) can be chosen arbitrary. Note that the l/e/^ 
divergences in soft diagrams do not always correspond to IR divergences in QCD, which 
is further evidence of their unphysical nature. This is the case ii D ^ B — C such as for 
the two loop graphs contributing to the running of Vk{v) 0. Finally, we see that setting 
^IR = ^uv in the soft amplitude in Eq. ( |A1| ) and running the ultrasoft modes from m to mv"^ 
with an anomalous dimension proportional to C — -B and running the soft modes from m 
to mv with an anomalous dimension proportional to A + B correctly performs the running 
between the scales. This is the method used here and in Refs. 

The correspondence in Eqs. ( [AID ^nd ( |A^ ) also provides a useful calculational tool: if 
the UV divergences A + C in the soft diagrams are known, and the combination B — C 
is determined from the UV divergences in the ultrasoft calculation, one arrives at {A + 
C) + {B — C) = A + B which is the combination needed to determine the soft anomalous 
dimension from A^. 



APPENDIX B: FUNCTIONS THAT APPEAR IN A£;NNLL 

The following functions of the principal and orbital quantum numbers {n,l,j) were de- 
rived in Refs. |23J23] and appear in our result in Eq. (HOl): 



^This addition might seem strange since it involves canceling an UV and IR pole. However, 
for a multiscale problem what is ultraviolet and what is infrared is always relative; if we 
label the e's by the corresponding scale then the cancellation is just l/e{mv) — l/e{mv) = 0. 
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No{n,l) 



N2{n,l) 



{Sl2)ljs 



+ / + 1) r 



^(n + / + l)-2 



nT{n - /) 



r(j + 2/ + 2) 



2r(n + / + l) f^^ r(j + l)(j + / + l 



n) 



r(j + 1) 



nr(n + / + 1) ^ 

+ 2r(n-/) ^.£^_^r(j + 2/ + 2)(j + / + l-n)3 

+ / + 1 



^(n + / + 1) + 7b] V ^'(^^ ^ ^ ^ ^ T2 

2r(n-/) r(2/ + 2 + j) 

"^r(r2 + / + i) f^, r(j + i)(j + / + i-n)2' 



f 2(/ + l) 
1 - 2/ 



-2/ 



2/ + 3 



: J =1-1 

■■3 = 1 

: j=/+l 



(Bl) 
(B2) 



(B3) 



(B4) 



Xijs = l{j{j + l)-l{l + l)-s{s + l) 



(B5) 
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